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Abstract 

Motivated by the excellent agreement between next-to-leading-order pQCD calculations 
and the inclusive jet spectra measured by CDF, we cross-check PYTHIA fragmentation func- 
tions. The convolution of the measured jet spectra with unmodified PYTHIA fragmentation 
functions results in reasonable agreement with the PYTHIA charged particle spectrum over 
the entire pr range of interest, while there is a sizable disagreement with the measured charged 
particle spectrum above pt — 30 GeV/c. In an attempt to understand the source of this dis- 
crepancy, we introduce a number of increasingly different toy-model fragmentation functions 
for the convolution. However, even the most extreme fragmentation functions result in an 
underestimate of the high-py CDF spectrum, which remains irreconcilable with the measured 
jet spectra. 



1 Introduction 



After an intense period of commissioning and a short period of collisions at 0.9 and 2.36 TeV 
in 2009 p], the continuous operation of the Large Hadron Collider (LHC) at multi-TeV collision 
energies looms just around the corner. Among the first measurements at the LHC will be the 
inclusive production of single charged particles (or hadrons), i.e. pp — > h + X, measured differ- 
entially in pseudorapidity (77) and transverse momentum {pt) [2]. While both distributions are 
the subject of Quantum Chromodynamics (QCD), the former is generally modeled phenomeno- 
logically due to the non-perturbative nature of \ow-px bulk production, while the latter, which 
involves hard processes, is generally described by the perturbative theory of QCD (pQCD). Hard 
production of high transverse momenta particles {px > 2 GeV/c) originates from the fragmen- 
tation of hard-scattered partons [3]. The fragmentation of hard-scattered partons into hadrons 
is described by the probability of finding a hadron carrying a specific fraction of the parton mo- 
mentum, known as the fragmentation function (FF). In hadronic collisions, a full description also 
requires knowledge of the distribution of the initial partons within the colliding hadrons, known 
as the parton distribution function (PDF). The measurement of the inclusive charged particle pt 
spectrum at large transverse momentum, therefore, measures in essence the convolution of three 
pieces: the hard-parton scattering cross section, the PDFs and the FFs. 

Especially at LHC energies, where a large fraction of the total cross section is comprised 
of the underlying QCD dynamics, a precise understanding of the QCD background rates is not 
only important for understanding Standard Model particle production (W^, Z, Higgs), but also 
for rare processes beyond the Standard Model [HE]. In addition, the inclusive charged particle 
Pt spectrum in pp collisions is an important reference for studying high-pj- particle suppression 
in the dense QCD medium produced in high energy nucleus-nucleus (AA) collisions [HIT]. The 
suppression (or enhancement) of high-pr particles is typically quantified by the ratio of charged 
particle pt spectra in AA collisions to those in pp collisions scaled by the number of binary nucleon- 
nucleon collisions, known as the nuclear modification factor Raa [Z]- At RHIC, the factor of 5 
suppression seen in Raa 13 El EH E] was an early indication of strong final-state medium effects 
on particle production. It is similarly expected to be one of the first measurements performed by 
the heavy ion programs at the LHC |12j . 

Experimentally, the inclusive charged particle pt spectra have been measured in pp and pp 
over a wide range of center-of-mass energies from 31 GeV to 1.96 TeV [121 El [121 ESI E]) and 
recently at 2.36 TeV at LHC [HICE]. While the measurements up to ^fs =1.8 TeV (and y/s =2.36 
TeV) are limited to px < 20 GeV/c, the latest CDF measurement [T7] at y/s =1.96 TeV, based on 
an integrated luminosity of 506 pb~ l , extended the pt reach up to about 140 GeV/c for the first 
time. However, the measurement shows that the high-p^ region cannot be described by the power- 
law modeling established in their earlier measurement at y / s=1.8 TeV |16j . The incompatibility 
between the new measurement and the power- law modeling grows from 50% at pt = 20 GeV/c 
up to a factor of 1000 at pr = 125 GeV/c. In order to fit the entire pt range, a more sophisticated 
parameterization was introduced [16] , namely another power law term was added to the previous 
fit. With the new parameterization, the normalized chi-square (x 2 /ndf) is reduced from 258/182 
to 80/223, albeit with a factor of 5-7 discrepancy remaining above ~ 90 GeV/c. The similarity 
between the spectra observed at both collision energies up to pr = 9 GeV/c (the range measured 
at 1.8 TeV) suggests that the incompatibility cannot be accounted for by the 9% increase in 
center-of-mass energy. 
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A disagreement of similar magnitude is observed when the measured high-p^ spectra are 
compared to leading-order (LO) and next-to-leading-order (NLO) pQCD calculations, as well as 
a simple extrapolation based on xt scaling [20, 21J. The validity of the factorization theorem 
is even questioned in [21], as a response to the huge discrepancy seen at high px - the regime 
where their NLO pQCD calculation should be most reliable. While the CDF paper [IT] offers no 
possible physics origin for the exceptionally large measured cross section at high p?, the sizable 
incompatibility with not only their former power-law modeling but also (N)LO pQCD calculations 
and xt scaling suggests that further study might be necessary. 

At the LHC, the nominal heavy ion collisions (PbPb) will take place at a center-of-mass 
energy of 5.5 TeV per nucleon pair, corresponding to the nominal pp collision energy of 14 TeV 
for the same magnetic rigidity Q For the first year, however, as the center-of-mass energy of 
pp collisions will be limited to 7 TeV [I], the corresponding PbPb center-of-mass energy will be 
limited to 2.76 TeV per nucleon pair. Since the Raa measurement requires a pp reference at the 
same collision energy as PbPb, the first-year heavy ion measurements at 2.76 TeV will rely on 
a combination of theory predictions and interpolations between lower energy measurements and 
those performed at 7 TeV. In this perspective, the CDF measurement is unique in two regards. 
First, the center-of-mass energy of 1.96 TeV is closest to that planned for the first-year heavy ion 
run (except for the measurement at 2.36 TeV which is limited in p? reach). Second, the reach to 
high pt far exceeds any previous measurements. Thus, understanding the observed discrepancy 
with the pQCD prediction is crucial. 

In contrast to the observed discrepancy in the charged hadron spectra, the CDF inclusive 
jet spectrum are in fact well described by NLO pQCD calculations (22J EH [Ml [25]. This is of 
particular interest, as high-p-p charged particles are understood in pQCD to be predominantly^] 
the fragmentation products of hard-scattered partons from the collision, i.e. "jets". In fact, 
NLO pQCD calculations carried out for charged particle spectra differ only from those for jet 
spectra by the addition of a parameterization of jet fragmentation. It is shown in [20] that the 
uncertainties related to different parameterizations of fragmentation functions (FF) and parton 
distribution functions (PDF) only amount to 10% and 25%, respectively. In this case, the only 
possible explanations for the exceptionally large measured cross section are either that the current 
modeling of fragmentation is dramatically incorrect (by a factor of 1000!), or there is a flaw in 
the measurement, or there is a breakdown in the QCD factorization theorem as suggested in |21| . 

In this paper, we attempt to reconcile the apparent discrepancy in the CDF inclusive charged 
particle spectrum by convoluting the CDF inclusive jet spectra measurement with a set of increas- 
ingly different fragmentation functions. We start with the PYTHIA fragmentation functions to 
see what the CDF jet spectra imply for the charged hadron spectra absent any surprises in the 
fragmentation. Then, we modify the fragmentation functions arbitrarily within the bounds of 
energy and momentum conservation to see if the measured hadron spectrum can be recreated. 
Finally, by using the hardest imaginable fragmentation function - each jet fragments into a single 
charged hadron - we rule out the possibility that unexpectedly hard fragmentation is responsible 

1 For the same magnetic rigidity of the LHC machine, the center-of-mass energy per nucleon pair in heavy ion 
collisions is just defined as the center-of-mass energy in pp scaled by the charge-to-mass ratio of the lead ion: 82/208. 

2 There is a contribution from the leptonic decays of weak gauge bosons, but this sub- leading processes is negligible 
in the inclusive spectrum. In PYTHIA it amounts to only 10% at most [20] ■ It is also possible for high-pr hadrons 
to originate in the absence of a jet in a higher-twist (HT) picture. However, as the production of hadrons in HT is 
power-law-suppressed in pr, that contribution should be negligible as well [26] ■ 
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for the exceptionally large measured cross section. 

The rest of the paper is organized as follows. In section 2, we introduce the technical 
details of the convolution method. We then describe the measured inclusive jet spectra and the 
PYTHIA inclusive jet spectra in section 3. In section 4, we argue that the convolution method 
can effectively reproduce the charged particle spectra within PYTHIA, i.e. by convoluting the 
PYTHIA jet spectrum with the PYTHIA fragmentation functions. We then show that the CDF 
jet spectra convoluted with the same PYTHIA fragmentation functions gives a similar result. In 
section 5, we compare to the CDF charged hadron spectrum the results of convoluting the CDF 
jet spectra with an arbitrary set of increasingly unrealistic fragmentation functions. The final 
section contains our conclusions and a further discussion on how this convolution method may be 
used to cross-check charged particle cross sections. 

2 Convolution Method 

In the QCD factorization scheme of hadron-hadron collisions [31 |27j, the invariant cross section 
for inclusive high-j»T hadron production is given by: 

d i a(AB -» CX) 
d 3 p c 

~ E /'Ik- C dx bq A a (x a -Q 2 ) q ^(x b ;Q 2 ) - DUz;Q 2 ) dHab ^ Cd) , (1) 
7T Jo io z at 

where the parton distribution function q A (x a ) describes the number density of constituents a 
within hadron A with longitudinal momentum fraction x a (in the range x a — > x a + dx a ). The 
fragmentation function B) c c (z) represents the probability that parton c hadronizes into C carrying 
a fraction z of the parton energy. Q 2 is the characteristic energy scale of the hard scattering. 
The LO cross section for the hard scattering of partons a and b at short distance is denoted by 
o. The summation is over all partons a, b, c, and d. The hadronization of parton d is implicit in 
the summation. A direct calculation of Eq. [T]is possible up to a certain order in a s provided that 
qa(xa) and D c c (z) are given. 

Similarly, the related cross section for inclusive jet production {AB — > Jet + X) can be 
given by the same equation, but with the term (l/z) D c c {z\ Q 2 ) replaced by 5(1 — z), since 
J2c Dc( z ) = <5(1 — z )- The convolution method is essentially equivalent to evaluating Eq. [ljbut 
based on measurements of the inclusive jet cross section and the fragmentation functions, provided 
that both are measured: 

Chad = PDF a / A (8) PDF b / B (g> a (hard parton scattering) (g> FF (2) 
= a jet (8) FF. 

In this case, it is clear that a knowledge of the inclusive jet cross section and the fragmentation 
functions associated by jet-pr are enough to reproduce the hadron spectra. The convolution of 
the jet cross section weighted in each jet-py bin by the associated fragmentation function can be 
cast into the following simple differential form: 
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where the summation is over all jet-pr bins, and the custom-built fragmentation function FFi is 
defined as: 



^ da had ' 

bbi{p T ,p Tjet ) = 



Act 



jet 



(4) 

PT,j e t=P^ jet 



The quantity FFi is just the transverse momentum differential cross section of charged particles 
per jet cross section at a certain jet-pr with finite bin size. 

The first term in Eq. [3] is well known not only from theoretical calculations but also from 
measurements of the inclusive jet cross section. However, the fragmentation functions are not 
known from any measurement in the exact form that they are needed (i.e. Eq. [4j). In particular, 
the convolution requires that FFi be measured in the same bins of jet-pr as the inclusive jet cross 
section measurement. 



3 Inclusive Jet Spectra (Data vs PYTHIA) 

The inclusive jet cross sections have been reported on numerous occasions by the CDF col- 
laboration, showing good agreement with NLO pQCD predictions for different jet algorithms^] 
[221 [221 [231 122] • ba particular, their latest measurement of the jet cross section [22], which used a 
mid-point cone algorithm on 1.13 fb^ 1 oipp collision data, agrees over a large range of jet-pr and 
rapidity not only with the NLO pQCD prediction within the respective experimental and the- 
oretical uncertainties, but also with the previous CDF measurements using different jet-finding 
algorithms. Good agreement with NLO pQCD predictions is also seen by the D0 measurement 
over a similar jet-pr range |28j . While a direct comparison of the inclusive jet measurements be- 
tween the two experiments is not possible due to different rapidity binning, systematic data-theory 
comparisons for the determination of LO, NLO, and NNLO pQCD PDFs show no significant in- 
consistency between the measurements |29j . In order to test whether the disagreement [17} [20] in 
the charged particle spectra between CDF and PYTHIA might originate from a difference in the 
respective jet cross sections, we compare the latest inclusive jet cross section measurement with 
that from PYTHIA. 

To maintain similar statistics over a large range of jet-pr, several QCD jets samples were 
generated in bins of the hard parton momentum transfer {pr = 0-15, 15-20, 20-30 GeV/c, etc.) 
using the D6T [3D] tune of PYTHIA 6.41 (or 6.42) [3T]. The different pr bins were then combined 
after properly weighting each sample by its corresponding cross section [^J Several different jet 
finding algorithms were used at MC level to test for possible algorithmic dependences of the jet 

3 kr [32] and midpoint cone [33] algorithms are used. 

4 (MSEL = 1) with low pr process (ISUB = 95) added for pr — > to avoid a divergent jet cross section. 
5 The cross section for the first pr bin was obtained by taking the difference of the minimum bias cross section 
with the sum of the cross sections for all other pr bins. 

Iterative cone, kr, and SIS cone were used with a cone radius R = J 'A</> 2 + Ay 2 = 0.7 
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cross section. The variation was observed to be smaller than the uncertainties in the measurement 
over the entire jet-pr range. 

Figure Qa) shows the inclusive jet cross sections from CDF data and PYTHIA, where the 
different rapidity (y) intervals have been scaled by arbitrary factors for clarity of comparison. 
In Fig. [ijb), where the different y bins have been plotted on the same scale, it is clear that 
the jet cross section decreases towards more forward rapidities. In the bottom of the figure, the 
ratio of the CDF jet cross section to the PYTHIA cross section is shown for each y bin. The 
inclusive jet cross section obtained from the PYTHIA D6T samples tends to underestimate the 
measured jet cross section. We note, however, that the default if-factor of cr NLO /a LO = 1 was 
used in the generation of our PYTHIA samples. Therefore, one may claim that the difference 
would be reduced by applying a greater-than-unity .ff-factor to account for the known difference 
between <r NLO and a LO [^] Other than the underestimation of the overall scale, PYTHIA describes 
the rapidity dependence of the jet spectra quite well for a wide range of jet-p-p. This study shows 
that the inclusive charged particle px spectra should not differ by more than a factor of 3 due to 
differences in the measured and PYTHIA jet cross sections (see lower panel of Fig. fTVb) ) . 



4 Convolution with PYTHIA fragmentation functions 

4.1 PYTHIA fragmentation functions 

The fragmentation functions defined in Eq. [4] are obtained from the same PYTHIA samples 
described in the previous section for the jet cross sections. Since these samples were generated in 
a series of pr bins, some care must be taken to ensure that the contributions from each pr sample 
are properly accounted for in both the numerator and the denominator of Eq. [4} Specifically, the 
fragmentation functions are determined via the following sum: 



dprdyhad \\yhad\<l-0) 



E 



A(J jet\\y Jet \ <y > jet prj<!( 
A (d N had/ d PTdyhad, 



N 3 In 3 
1 * event! u e 



event 




1 Y event I u ei 



event 



(5) 



where N d event and cr 3 event are the number of events and the cross section for the j 



:ifl pr bin, re- 
spectively. Each fragmentation function is evaluated at prjet = Pr j e t w hh a width of Ap l T - et for 



the i jet-j?T bin. Nhad an d ^jet are the number of charged particles and the number of jets in 
\yhad\ < 1-0 and \yj e t\ < v'j e ti respectively. (Note that yhad and yj e t are not necessarily the same.) 

The PYTHIA fragmentation functions FFi(pT,p l T Jet ) for \yhad\ < 1-0 and \yj e t\ < 1-0 are 
shown in Fig.[2|a) for px < 140 GeV/c and bins oip % T j et corresponding to the CDF jet cross-section 
measurement. The hardening of the fragmentation function with increasing jet-pr is apparent. 
The same fragmentation functions are shown in Fig. [2^b), where each has been weighted by 
the associated jet cross section in order to compare their relative contributions to the charged 



7 For example, the A'-factor for the inclusive charged particle spectra at Tevatron energy (^fs =1. 
phenomenologically to be slightly above unity (1.28) [34] 



TeV) is found 
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Figure 1: (a) The inclusive jet cross sections in different rapidity intervals. CDF values (filled 
markers) measured at hadron level using the midpoint algorithm (cone radius=0.7); error bar is 
statistical only. PYTHIA generated with D6T tune (solid lines) using an iterative cone algorithm 
(radius=0.7). An arbitrary factor of 10 3 separates the different rapidity intervals as in the CDF 
paper for clarity, (b) The same distributions as in (a), but with all rapidity bins on the same 
scale. The ratios of the PYTHIA jet cross sections to the measured values are plotted in the lower 
panel with systematic uncertainties drawn for the measurement in 0.1< \y\ <0.7 to illustrate the 
size of the uncertainties involved in the measurement. 
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particle spectra. Due to the steeply falling nature of the jet cross section, the vast majority of 
charged particles in the region of interest (30 GeV/c < px < 140 GeV/c) are from jets with 
Pt < 400 GeV/c. The contribution from higher px jets is less than 0.1%. 




Figure 2: (a) PYTHIA D6T fragmentation function in each jet-pr bin. (b) The same fragmenta- 
tion functions weighted by the associated jet cross section measured by CDF. da CDP corresponds 
to da j e t I dpxdy in Eq. [3j 

Although a direct verification of the PYTHIA fragmentation functions in Fig. [2] is not 
possible in the absence of the same measurement on 1.96 TeV data, reasonable (and sometime 
very good) agreement |^] has been seen between PYTHIA and a variety of fragmentation-related 
measurements [35l [361 [371 EH]- For example, detailed CDF studies of inclusive jet shapes at 
-y/s=1.96 TeV [38] show that the integrated and differential jet shapes are well described by 
PYTHIA Tune A for jet-pr up to 380 GeV/c. Furthermore, the fact that PYTHIA and NLO 
charged particle spectra are in reasonable agreement [20J, where the latter uses fragmentation 
functions based on global data fits (AKK [39J, DSS 00], HKNS [41]), implies that PYTHIA 
fragmentation should be comparable to these global fits. In this case, the fragmentation functions 
obtained from PYTHIA should be a reliable proxy for the global understanding of fragmentation 
from data. Regardless, this paper investigates discrepancies significantly larger than 50%, so 
the detailed matching of the PYTHIA fragmentation model to measurement is not our primary 
concern. Moreover, model-independent fragmentation functions will be introduced later in Section 
5. 

8 There is some indication that the fragmentation properties of quark jets are rather poorly described by PYTHIA 
when quark and gluon jet are investigated separately '35 . 



8 



4.2 Convolution of PYTHIA jet spectrum with PYTHIA fragmentation functions 



To test the effectiveness of the convolution technique put forth in Section [2j we first attempt to 
retrieve the known PYTHIA charged particle pr spectrum from a convolution of the PYTHIA jet 
spectra with the PYTHIA fragmentation functions calculated in Eq. [5] and plotted in Fig. [2j We 
perform the convolution separately for each of the five jet rapidity intervals (see Fig. [TJ, in order 
to quantify their relative contributions to the single particle spectra. The results of this test are 
shown in Fig. [3j where the PYTHIA charged particle spectrum (solid line) is compared to the 
output of the convolution for charged particles with \n\ < 1.0. The relative contribution of jets in 
different y ranges to \n\ < 1.0 charged particles can be seen in Fig.^a), where the rapidity ranges 
are indicated by the same symbols as in Fig. [I] As expected from the measured jet shapes^] [38] . 
the large majority of particles fragment from jets within \y\ < 1.1. In Fig. [3^b), the contributions 
from all rapidity intervals are summed and compared to the true charged particle spectrum. In 
the bottom of that figure, the ratio of the convoluted spectra to the true charged particle spectra 
is shown. Except at low pt (below a few GeV/c) where non-perturbative particle production 
and jet reconstruction inefficiency become relevant, the convolution method reproduces the true 
spectra almost exactly (well within 0.01% for pr > 6 GeV/c). 



4.3 Convolution of measured jet spectrum with PYTHIA fragmentation functions 

With the robustness of the convolution technique verified on PYTHIA, the next step is to introduce 
the CDF measurement of the jet spectrum in the convolution. Given the level of agreement already 
demonstrated in Section [3] between PYTHIA and the measured jet spectra, one should expect the 
convolution based on the measured jet spectra to result in quite similar charged particle spectra 
as already seen for PYTHIA in Fig. [3j However, there is an additional complication involved in 
using the measured jet spectra in the convolution instead of the generated PYTHIA events. 

Unlike for PYTHIA where the jet cross section is available down to very low pr, the jet 
cross section measured by CDF has only been published between 62 and 700 GeV/c. The ef- 
fect on the charged particle spectrum from excluding the fragmentation products of jets below 
62 GeV/c is shown for PYTHIA in Fig. |3^b) represented by empty squares. These low-pr jets 
contribute significantly to the charged particle spectrum up to around 40 GeV/c. As expected, 
the contribution completely vanishes at pr = 62 GeV/c, since a jet cannot fragment into a more 
energetic charged particle. Since the disagreement with (N)LO pQCD calculations is most promi- 
nent above pt > 50 GeV/c, a detailed understanding of the low-pr contributions is not central to 
this investigation. Hereafter, whenever the measured jet cross section is used in the convolution, 
we use the PYTHIA value for pr < 62 GeV/c. 

In Fig. |4](a), the charged particle pr spectra are shown for the convolutions based on the 
measured jet cross sections, again using the same symbol conventions for the rapidity ranges as 
in Fig. [I] In Fig. [4](b) , the contributions from all rapidity intervals are summed as it is done 
in [3^b). To see the variation of the obtained spectrum due to the uncertainties in the CDF 
jet measurement, a conservative choice of 82%, which is the largest systematic uncertainty from 
the measurements in the first three rapidity intervals for }et-pr below 457 GeV is applied to 
the obtained charged particle spectrum, which is shown as a grey band. Also, a grey band is 

9 From measured jet shapes, we know that the majority of the energy in a jet is concentrated around jet axis. 
So, the rapidity of the leading tracks should correlate quite closely with the rapidity of the jet. 
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Figure 3: (a) PYTHIA charged particle cross section obtained from the convolution of the 
PYTHIA inclusive jet cross sections in different jet y ranges with PYTHIA fragmentation func- 
tions (filled markers) compared to the "true" charged particle cross section from PYTHIA (solid 
line), (b) The convoluted spectra after summing the contributions from all jet y ranges (filled 
circles), and the same after excluding the contribution from jets with pr < 62 GeV/c (open 
squares). In both figures, the CDF measured cross sections are also shown (empty circles) for 
comparison. 
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drawn for the measured charged particle spectrum to indicate the size of the uncertainties in 



the measurement 10 The spectra from the PYTHIA-only convolution is drawn as black lines for 
comparison in bothQa) and[4]^b). In the bottom of the figure, the ratio of the resulting spectrum 
from the convolution of PYTHIA jet spectra to that of CDF jet spectra is shown. As expected 
from the agreement between the CDF measured jet spectra and PYTHIA (better than ±50% 
below 400 GeV), the convoluted charged particle spectra are in similar agreement. 




Figure 4: (a) Charged particle px differential cross section obtained from the convolution of 
the CDF jet cross sections measured in five rapidity intervals with the PYTHIA fragmentation 
functions (filled markers). Same cross sections obtained from the PYTHIA-only convolution are 
also shown in comparison (black line), (b) The convoluted spectra after summing the contributions 
from all jet y ranges (filled circle), and the same from the PYTHIA-only convolution (black line). 
The ratio of the two resulting charged particle spectra (from convolution with either the CDF or 
PYTHIA jet spectra) is shown in the lower panel. 



5 Convolution with toy-model fragmentation functions 

Knowing that the charged particle spectra from the convolution of the measured jet spectra with 
PYTHIA fragmentation functions still vastly undershoot the measured spectrum, it is interesting 
to know how sensitive the charged particle spectra are to arbitrary changes in the fragmentation 

10 Statistical uncertainty only, although the statistical uncertainties are comparable to the total uncertainties for 
p T > 50 GeV/c |T7| 
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function; or rather, how large a change in the fragmentation functions would be required to 
reproduce the measured spectrum? We address this question by convoluting the measured jet 
spectra with the fragmentation functions obtained based on the following two toy fragmentation 
models: 

1. Harder fragmentation - modified shapes 

2. Hardest possible fragmentation - one charged particle per jet 

The modified fragmentation functions are obtained in the same manner as for the default PYTHIA 
(see Eq.[3]). 

5.1 Harder fragmentation 

We make the default PYTHIA fragmentation functions harder by modifying the shape of the 
fragmentation functions. Of course, there are arbitrarily many imaginable functional forms one 
could use for "hardening" the fragmentation functions. Here we have chosen to harden the 
fragmentation functions by bending the pr-shape with a simple power-law functional weight: 

f(p T ) = 1 + c(p T - pT n T for PT > P T n , (6) 

The shape of the modification is governed by the constant c and the exponent re, while pji m 
determines the px value where the modification begins. To maximize the effect in this study, we 
bend the shape of the fragmentation function to the extent that the energy sum of the fragmented 
charged particles equals the energy of the corresponding jet. In this fashion, the fragmentation 
function can be maximally hardened for any choice of two parameters: c and p™ m . We try four 
different combinations of c and p^p m , which are presented in Fig. [Ej The values of c (0.005 and 
0.0005) were selected to ensure the bending is not unrealistically abrupt. Since the spectrum 
from the convolution already describes the CDF measurement rather well up to 20 GeV/c, we 
pick p™ m accordingly (10 and 20 GeV/c) to maintain this agreement. With these choices of c 
and p™ m , the exponent n that fulfills the maximum hardening ranges from 1 to 4 for the various 
jet-pT bins. The modified fragmentation functions from jets within a range of 0.1 < \y\ < 0.7 
are shown in Fig. [5] for the different values of c and p™ in . Here, we only show the fragmentation 
functions for jets with pt up to 400 GeV/c, as the fragmentation functions with higher jet pt have 
a negligible contribution to the inclusive spectrum (see Fig. [2^b)). When compared to Fig. [2^a), 
the hardening of the fragmentation functions is immediately apparent. Despite an attempt to 
maintain smooth functions, the fragmentation functions in the first few jet-py bins are if anything 
unrealistically hard, as they contain regions with positive slopes. 

For each of the four scenarios, the convolution of the hardened fragmentation functions 
with the CDF jet measurement is performed for the three dominant rapidity intervals (|y| <0.1, 
0.1< \y\ <0.7, 0.7< \y\ <1.1), which are then summed. In Fig[6^a), the resulting charged particle 
spectra are compared to the CDF measured charged particle spectrum. The harder fragmentation 
functions are reflected in a hardening of the inclusive spectrum. The disagreement between the 
measurement and the resulting spectra is less pronounced between 30 and 60 GeV/c. However, 
a disagreement of up to two orders of magnitude remains at the highest measured transverse 
momentum. 
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Figure 5: (a-d) PYTHIA (with D6T) fragmentation function in each jet-pr bin modified by the 
power-law function (see Eq. [gJ) for different choices of the "hardening parameters", c and p™™. 
The same color and line style conventions are used from Fig. [2} 
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5.2 Hardest possible fragmentation 

Finally, we assume the extreme case of hard fragmentation, for which each hard-scattered parton 
fragments into a single charged particle. This is the hardest possible fragmentation [z = 1) that 
still conserves energy (though not necessarily other conserved quantities). Unlike the previous 
fragmentation functions, which are still model-dependent to some extent, the z = 1 scenario is 
model-independent . 

In Fig [6^b), the result of applying the z = 1 fragmentation to the measured jet spectrum is 
shown. Again, the contributions from the three dominant rapidity intervals are summed to obtain 
the inclusive charged particle spectrum for \rj\ <1. Both the shape as well as the overall magnitude 
are dramatically changed, enhancing the cross section by a few orders of magnitude. In the same 
figure, the jet spectra for the three y intervals (\y\ <0.1, 0.1< \y\ <0.7, 0.7< \y\ <1.1) are averaged 
and plotted in the form of an invariant yield for comparison. However, the charged particle 
spectrum measured by CDF at the highest pt actually exceeds their measured jet spectra with 
z = 1 fragmentation. Furthermore, the level of underestimation is beyond the level of systematic 
uncertainties shown in the CDF jet measurement, which range from 10 to 80% depending on 
jet-pr M- 

6 Summary and discussion 

The latest CDF measurement of the inclusive charged particle px spectrum up to very high px is 
cross-checked against the convolution of the measured inclusive jet spectra with a set of different 
fragmentation functions. The inclusive charged particle spectrum obtained from the PYTHIA 
fragmentation functions convoluted with the measured jet cross section fails to reproduce the 
measured spectrum, despite matching (N)LO predictions reasonably well. PYTHIA fragmentation 
functions modified by a simple toy model result in quite different spectra shapes. However, despite 
an improved agreement at intermediate px, the convolutions still underestimate the measured 
charged particle spectrum at high pr- Finally, we show that even the most extreme case of each 
jet fragmenting into a single charged particle fails to reconcile the measured jet and charged 
particle spectrum. Based on these studies, we rule out the possibility that the disagreement of 
the NLO pQCD calculations with the measured charged particle spectrum (but not the measured 
jet spectra) is due to an unexpectedly hard fragmentation of jets. Thus, we conclude that the 
CDF charged particle spectrum cannot be reconciled with the present understanding of factorized 

p qcd[5 

This work was originally motivated for the purpose of cross-checking the CDF charged par- 
ticle pt spectrum measurement, but this convolution technique can be extended for use in future 
measurements without relying on input from PYTHIA at all. Any of the following three measure- 
ments can be independently obtained from the (de)convolution of the other two: the inclusive jet 
spectrum, the fragmentation functions of inclusive jets, and the inclusive charged particle spec- 
trum. We expect that such a cross-check of future single particle px spectra measurements against 
the available calorimetric information will be very helpful for understanding the tracking-related 
systematic uncertainties in the high-p^ regime. 

11 This doesn't rule out the possibility discussed in [35] that the factorization breaks down, since this investigation 
with the convolution method relies upon the validity of the factorization theorem. 
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Figure 6: (a) Charged particle pt differential cross sections obtained from the convolution of the 
measured jet cross sections with PYTHIA fragmentation functions modified according to Eq. [6j 
Different symbols show the spectra obtained with the modified fragmentation functions based on 
different sets of hardening parameters, (b) The same cross sections obtained from the convolution 
of measured jet cross sections (|y| <1.1) with the hardest possible fragmentation (filled circles) 
compared to PYTHIA (solid line). The CDF combined jet spectra averaged over |y| <1.1 (filled 
triangle) and the CDF charged particle spectrum (open circle) are plotted for comparison. 
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Note added: Another analysis reaching similar conclusions [32] was submitted for publication 
during the finalization of this paper. 
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